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Abstract: We study the center symmetry of SU(A^) gauge theories with fermions in the 
two-index representations, by computing the effective potential of the Polyakov loop in the 
large-mass expansion on the lattice. In the large- limit and at non-zero temperature, we 
find that the center symmetry is for fermions in the adjoint representation and just Z2 
for fermions in the (anti) symmetric representation. We discuss the fact that our results do 
not contradict the orientifold planar equivalence, which relates a common sector defined by 
the bosonic gauge-invariant C-even states of theories with fermions in different two-index 
representations. Our results complement the work of Armoni et al. (2007), who showed 
how at zero temperature a center symmetry is dynamically recovered also for fermions 
in the (anti)symmetric representation, by considering the theories at finite temperature. 
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1. Introduction 

Confinement of colour is related to center symmetry: where center symmetry is present and 
not broken, it implies confinement. Its spontaneous breaking leads to the deconfinement 
transition, and the critical behaviour at the phase transition is described by a three- 
dimensional effective theory which encodes the relevant pattern of symmetry breaking [l| . 

On the other hand, orientifold planar equivalence is the equivalence in the large-A^ 
limit of the SU(A^) gauge theory with a Dirac fermion in the symmetric/antisymmetric 
two-index representations (OrientiQCD) and the SU(A^) gauge theory with a Majorana 
fermion in the adjoint representation (AdjQCD) in a common sector ||2|, |3|, The validity 
of this equivalence was discussed in detail in Refs. ^. A proof of the equivalence on 
the lattice in the strong-coupling and large-mass phase can be found in Ref. [0] (a general 
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setup for descussing planar equivalences between theories with two-index representations 
in the strong-coupling and large-mass phase was also presented in |pl H ) . 



A potential inconsistency was pointed out in Ref. |1C]: the two theories involved in the 
equivalence have different symmetries and in particular different center symmetries. The 
action of the adjoint theory is invariant under the full center Zj^ of the gauge group at 
every value of N. In the orientifold theories the center symmetry is explicitly broken by 
the matter to a Z2 N is even and to nothing if N is odd. The different symmetry pattern 
in the two theories suggests a different dynamics at the deconfinement phase transition. 
It is therefore interesting to investigate to which extent the orientifold planar equivalence 
can describe correctly the dynamics at the deconfinement phase transition. 



A solution to this puzzle in the case of zero temperature was presented in Ref. |11|: 
in the orientifold theories the action can be separated in a part that is invariant under 
Zjv and a part that is not. In the large- limit, the non-invariant part decouples from 
the expectation values of the Wilson loops. The Z^v symmetry is therefore dynamically 
recovered and the apparent contradiction is resolved. One of the physical consequences is 
the stability of the fc-strings. 

In this work we consider the orientifold theories at finite temperature, both in the 
confining and deconfining phases. The main difference with the zero temperature case is 
that now the theory contains also loops wrapping around the thermal direction. We show 
by explicit computations that the non-invariant part of the action decouples only from the 
expectation values of the loops with zero winding number around the thermal direction. In 
the general case of wrapped loops, the non-invariant part of the action does not decouple, 
thus the Z^r symmetry is not restored in all the sectors of the theory. In particular, we 
will show that matter in the orientifold theories explicitly breaks the center symmetry to 
the Z2 subgroup in the large-A limit at finite temperature. 

In the deconfined phase, this result is well known. It can be obtained in the high- 
temperature regime, by analysing the one-loop effective action for the Polyakov loop. The 
invariance under Z2 but not Zjv implies that in the deconfined phase (where the Z2 is 
spontaneously broken), only two degenerate vacua exist. 

In the large- A limit, the gauge theories discretised on the lattice exhibit two confined 
phases as the 't Hooft coupling, A = 5^ A, is varied: an unphysical strong-coupling phase 



and a physical weak-coupling one [l^. In both these phases, we investigate the sym- 
metry content by computing the functional generator (W) of the connected expectation 
values of the Polyakov loop in the framework of the large-mass expansion. The generator 
is a legitimate tool to study the center symmetry, since it encodes the same amount 
of information as the effective potential of the Polyakov loop. For finite A^, W can be 
expanded in powers of 1/A: 

W(A,Ar) = ^a„(Ar)A-" . (1.1) 

n 

At strong coupling the large-A limit of W is obtained by taking the limit of each coefficient 
an{N) in the above expansion. Properties of Polyakov loops correlators can be inferred 
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from W. In particular we will show that: 



hm ((tr 0)2), / , (1.2) 

N-^oo 

in the orientifold theories, where (tr 0,) is the Polyakov loop; the above result shows clearly 
that the full center is not a symmetry. We may add that it is explicitly broken (and not 
spontaneously), since the vacuum is not degenerate in the strong-coupling regime. Under 
some mild assumptions, this result remains valid in the physical (weak-coupling) confined 
phase. 

Let us now anticipate some comments on the original question. How can the orientifold 
planar equivalence be valid if the two theories that should be equivalent have got different 
symmetries? It was already pointed out that the equivalence holds only in a common 
neutral sector of the various theories. The neutral sector was identified in Ref. [^] to be 
the set of all the single-trace gauge-invariant observables that are also invariant under 
charge-conjugation (C) symmetry. The Polyakov loop does not belong to the neutral 
sector; indeed under C-symmetry (tr fi) — s- (tr Q)^ , therefore the real part is C-even and 
belongs to neutral sector, while the imaginary part is C-odd (it belongs to the twisted 

I'Kik 

sector). Under the action of an element e'^ of Z^v, the neutral and twisted sectors are in 
general mixed: 

Re [tr ^ cos (^^^ Re [tr Vl] - sin (^^^ Im [tr SI] , (1.3a) 

/ 2TTk\ / 27rk\ 

Im [tr n] sin ( — J Re [tr J^] + cos ( — J Im [tr n] . (1.3b) 



The adjoint theory has a Zf^f symmetry, but the neutral and twisted sectors are mixed 
by the center symmetry. Only the Z2 subgroup (which corresponds to e ^ = — 1 in 
Eqs. ( |1.3| )) maps the neutral sector into itself. Therefore the orientifold planar equivalence 
requires just a Z2 symmetry in the neutral sector of the orientifold theories. The fact that 
Z2 (and not Z^) is the center symmetry for the full orientifold theories does not follow 
from the equivalence and is not in contradiction with it. 

The paper is organised as follows. In Sect. ^ we give a brief review of center symmetry. 
In Sect. 1^, and ^ both the strong-coupling and weak-coupling confining phases on 
the lattice are analysed. Although the general setting is valid for both phases, explicit 
analytical computations can be performed only in the strong-coupling phase. In particular, 
the effective potential in the large-mass phase is computed. Even though this phase does 
not describe continuum physics it is interesting because it provides a testing bed where the 
planar equivalence holds and computations can be performed explicitly. 

The possibility of extrapolating some of the results from the strong-coupling regime 
to the continuum limit is discussed in Sect. ^, together with the necessary underlying 
assumptions. A brief analysis of the center symmetry in the high-temperature regime is 
presented in Sect. ^. 

Finally in Sect. 0, we introduce the effective potential in the Hamiltonian formalism and 
investigate the center symmetry by means of the formalism of the coherent states [^, ^. 
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2. Center symmetry 



Whenever a gauge group has non-trivial center, the center symmetry plays a central role in 
understanding the dynamics at the deconfinement phase transition. We briefly summarise 
in this section some relevant properties of gauge theories under transformations that belong 
to the center of the gauge group. The same notation is used in the rest of the paper. For 
a gauge theory discretised on the lattice, a transformation belonging to the center of the 
gauge group acts on the link variables according to: 

Uo{xo = 0, x) ^ uUo{xo = 0, x) , (2.1) 

where xq is the coordinate along the time direction, and u is an element of the center. 
All the other link variables are left unchanged by the center transformation. The center 
transformation counts the winding number of Wilson loops around the time direction. If 
the closed path F wraps w times around the time direction, the transformation rule is: 

tr [/(F) ^ -u^'tr U{r) , (2.2) 

where U (F) is the path-ordered product of the link variables along F. Since the gauge 
action is given by a sum of plaquettes, which have zero winding number, it is invariant 
under the action of the center. 

A generic representation R of the gauge group induces a representation of the center 
subgroup, that can be labelled by an integer Nr (with < Nr < TV), the AT-ality of the 
representation: 

R[u] = ti^^ , R[uU] = u^^'RlU] . (2.3) 

Clearly, if matter in a representation with zero A^-ality is present, the theory is still in- 
variant under the action of the center. Otherwise, the theory is invariant only under some 
subgroup. The center symmetry of the theories investigated in this work are summarised 
in the following table. 



gauge group 


matter 


center symmetry 


S0(2iV) 


none/adjoint 




SU(iV) 


none/adjoint 


Zn 


SU(7V) 


fundamental 




S\J{N) with N even 


AS/S 




SU(iV) with N odd 


AS/S 




U(iV) 


none/ adjoint 


U(l) 


U(iV) 


fundamental 




U(iV) 


AS/S 


Z2 



3. Potentials for the Polyakov loop 

In this section, we introduce effective potentials that describe the dynamics of Polyakov 
loops for the gauge theory discretised on the lattice. We shall first discuss the case of a pure 
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gauge theory, while fcrmions in arbitrary representations arc introduced in the next section 
within the large-mass expansion framework. In what follows, it is useful to distinguish the 
temporal coordinate xq from the spatial ones, i.e. x = (a;o,x), and the temporal link 
variable Uo{xo,yi) from the spatial ones ?7fe(xo,x). 

The parallel transport CI around the temporal dimension and the Polyakov loop P are 
defined as usual by: 

Q(x) = C/o(0,x)i7o(l,x)---C/o(r-l,x) , (3.1) 
P(x) = tr n(x) , (3.2) 

where T is the extent of the lattice in the temporal direction in units of the lattice spacing. 
Gauge invariance is used to fix the temporal gauge, where we have: 

f 1 if xo = 0, . . . , T - 2 , , 

C/o xo,x = <^ _ • 3.3 

[n{:>c) ifxo = r-i 

After such a gauge fixing, the gauge action 5(J1, Uk) is invariant under time-independent 
gauge transformations. The Fadeev-Popov determinant being trivial, the measure of the 
functional integral is simply: 

e^p{-S{n,Uk)} H dUk{xo,^)l[dn{^) . (3.4) 

Let us now define several functional of the field ri(x); they all encode information 
about the dynamics of the Polyakov loops, and provide useful information on the structure 
of the deconfinement phase transition. 

The probability distribution of the parallel transport Jl(x) is obtained by integrating 
out all the link variables in the spatial directions: 

e-5o{n) ^ I y e-5(n,[/.) "Q dUk{xo,^) . (3.5) 

The functional Sq{^) defines an effective action for the Polyakov loop, and the vacuum 
expectation value of a generic functional /(fi) can be computed as: 

{fm = I / /(f^)e-^-(^) Hdni^) . (3.6) 

X 

The functional S'n(i^) induces a natural definition for the probability distribution of 
the Polyakov loop as the expectation value of the delta function: 

^-N^V{P) ^ (JJ<5(p(x) - tr 0(x))) = 

X 

= i y dn{x) Yl 6 (P(x) - tr 0(x)) e-^"(^) . (3.7) 
The function V{P) is sometimes referred to as the effective potential for the Polyakov loop. 
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In this work, we shall concentrate instead on the quantum action for the Polyakov 
loop, which is defined as the Legendre transform of the free energy of the system. First we 
construct the generator of connected n-point functions by coupling the Polyakov loop to 
a complex external source: 



g-7v2iy(a,a) ^ J rfj](x)e-^f^('^) exp |-iV [a(x)P(x) + a(x)pt(x)] | 
A Taylor expansion of the exponential in Eq. (^^) yields: 

oo ^ 

w{a, a) = 1 - ^ (-If — • • • «(yi) • • • «(y<i 



(3.8) 



p,q=0 



p\ql 



Xl,...,Xp 



(P(Xl)...P(Xp)Pt(y,)...pt(y^))^ 



where the (properly normalised) connected expectation values are defined as: 

<P(Xl)...P(x,)Pt(y^)...pt(y^)) _ 



(3.9) 



5a{-xi) ■ ■ ■ Sa{-Xp)6a{yi) ■ ■ ■ Sa{yq) 



W{a, a] 



(3.10) 



a=a=0 



Finally, the Legendre transform of the generator of connected Polyakov correlators yields 
the generator of IPI diagrams: 



z(x) 



-W{a, a) 



(P(x) exp{-(aP + aPt)}; 



a(x)"''^'~^' (exp{-(aP + aPt)}) 

r(z, z) = W{a, a) — {az + az) , 



(3.11a) 
(3.11b) 



where we have omitted the dependence on the spatial coordinates in order to simplify the 
notation. The functions z{a,a), z{a,a) yield respectively the expectation values of the 
Polyakov loop and its complex conjugate in the presence of the external sources a, a. We 
shall refer to the function r(2;, z) as the quantum action for the Polyakov loop. In general, 
it is different from the function V{z); however it can be shown that the two functions 
coincide for space-independent Polyakov loops in the limit of infinite spatial volume^. 

The center is a symmetry group of the quantized theory if and only if the effective 
potential T{z,z) is invariant^. Because of the Legendre transform, the effective potential 



^The provided definition of the quantum action and the decomposition of VK(q, a) in connected expecta- 
tion values are correct as long as the vacuum is non-degenerate. Otherwise the relationship between (a, a) 
and {z, z) is not invertible and the definition of the quantum action is more involved; while a sum over all 
the vacua is required in the Taylor expansion of W{a,a). The non-degeneracy condition is satisfied in the 
strong-coupling phase, in which we perform explicit computation. 

^In principle, one should consider the effective potential of a generic loop arbitrarily wrapping along the 
time direction. This can be done using the same formalism developed here. Ifowever we shall assume that 
the Polyakov loop completely characterises the center symmetry. 
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is not easily computed. However, since z{ua,ua) = uz{a,a), we can readily prove that, 
for any element of the center u: 

T{uz, uz) = W{ua, ua) — az{a, a) + az{a, a) , (3-12) 

and the following implications hold: 

center is a symmetry 44> T{uz, uz) = T{z, z) 44> W{ua, ua) = W{a, a) . (3.13) 

Therefore, the functional W{a,a) encodes all the information that we need about the 
center symmetry. The terms of the expansion in Eq. (|3.9D are classified in Tab. || according 
to their symmetry properties. 

Term in W Symmetry 
p = q mod 2 Z2 
p = q mod N 

p = q U(l) 

Table 1: Symmetry properties of the Polyakov connected correlators {{P)p{P'')'^)c. We say that 
G is a symmetry group for the Polyakov connected correlator if the latter is invariant under the 
transformation ri(x) ur2(x), with u G G. Note that G is not necessarily a subgroup of the center. 



The connected correlators can be computed analytically for the pure gauge theory in 
the strong coupling regime, where a topological classification of diagrams can be setup 
which closely follows the one obtained in the continuum by 't Hooft It]. The details 



relevant for our computations are summarised in App. The large-A^ behaviours of the 
correlators are: 

{{P)nP^)'^)c,YM 

^^2=^3^ =0(1), lfp = 9, 

{{P)nP^Y)c,YM , 

— C — , \i p = q mod A* , p q , 



0, ifp/g mod AT. (3.14) 



4. The connected graphs in the large— mass expansion 

Let us now discuss the introduction of fermions in the theory. Since we are interested in 
the orientifold planar equiavalence, we shall consider fermions in arbitrary representations. 
The fermionic effective action can be written as a sum of Wilson loops in the large-mass 



phase [17|: 



SfiUk, ^) = -Yl c(^) PPiu;)] , (4.1) 



where lo indicates a generic closed path on the lattice of length L{lj), U{uj) is the parallel 
transport along u), and R\U\ is the matrix representing U in the representation R. The 
coefficients c(a;) depend on the spin structure of the chosen discretisation of the Dirac 
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operator and are independent of the colour structure. They are of order of m~^^'^\ where 
m is the bare fermion mass. It will be convenient to split the set C of all the closed paths 
into the union of the sets C{w) of closed paths with winding number w around the thermal 
dimension. 

The connected expectation value of Polyakov loops can be formally thought as a func- 
tion of the coefficients c(u;). In the large-mass phase, it can be expanded as a power series 
around c(cj) = 0: 



n=0 u)i...LOu ^ ^' ^ 



Sf=0 



(4.2) 



Each derivative with respect to c{uj) corresponds to the insertion of a vertex tr R[U{uj)] 
in the expectation value. Since the derivative must be computed at Sf = 0, each term in 
the expansion above can be written as combinations of connected expectation values in the 
pure Yang-Mills, as shown in some explicit examples later on. Therefore the asymptotic 
behaviour of the connected expectation values in the theory with dynamical fermions can 
be reconstructed from the asymptotic behaviour of the connected expectation values in the 
pure Yang-Mills. 

Let wi, . . . be some closed paths and wi, . . .Wnhe their respective winding numbers 
around the thermal direction. Consider the generic normalised connected expectation value 
in the pure Yang-Mills theory: 

(Hjtr U{uJi))c,YM . 

This quantity is always finite in the large- limit (it could be also and in this case we say 
that the leading order is not saturated). This fact is proved in the strong-coupling, and in 
the perturbative expansion, and it is consistent with results from numerical simulations. 
Moreover, it is equivalent to the N'^ scaling of the free energy of the system, when external 
sources of order N are coupled to the operators tr U{u>i) in analogy to the definition (^), 
and consequently it is equivalent to the A^'^ scaling of the quantum action for the expectation 
values of the operators tr U{uji)/N. The following cases can arise: 

• Wi mod N. The expression in Eq. (^^) transforms non-trivially under center 
transformations. Since the center symmetry is not broken in the confined phase, the 
expression is exactly zero for every value of N. 



^iWi = mod N but X^jtfi / (see Figs, y). The expression in Eq. (O) is 
invariant under Z^r but not U(l). Using the strong-coupling expansion (see App. ^ 
it is possible to prove that the correlator vanishes in the large- A'^ limit. This result 
can be extended beyond the strong-coupling phase, under the commonly accepted 
assumption that the U(A^) and SU(A^) Yang-Mills theories differ in the large- A^ limit 
only for subleading contributions. If the gauge group were U(A^), the expression in 
Eq. ( f4.3D would be exactly zero for every value of A^. Since we are interested in the 
SU(A^) gauge group, it must be zero in the large-A^ limit. 
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The loops with winding number different from zero appear in pairs with opposite 
winding number, as shown in Fig. [l|. In particular Wi = 0. From the strong- 
coupling expansion (see App. ^) we get that the expression in Eq. ( |4.3| ) in this case 
is different from zero in the large- limit. We will assume that this result extends 
beyond the strong-coupling phase. 



Wi = but some of the loops with winding number w cannot be paired with 
a loop with winding number —w. Although the expression in Eq. (4.3) is invariant 
under the center (both Zjv and U(l)), it vanishes in the strong-coupling phase (see 
App. ^) because of the topology of the torus. We cannot say if this result extends 
beyond the strong-coupling expansion, but anyway we will not use it in this work. 







Figure 1: Representations of some connected expectation values of products of loops with net 
winding number equal to zero, in pure Yang-Mills. The first one is (tr ft tr fl^)c,YM'i the Polyakov 
loops are connected by an oriented surface that wraps around the thermal dimension. In the strong- 
coupling expansion, the surface is tiled by plaquettes coming from the Wilson action; products of 
two (or more) group elements associated to each link of the lattice are integrated with respect 
to the Haar measure. This graph yields a contribution 0{N'~'). The second graph represents 
((tr tr ri^)^)c,YM; the two tubes are glued together through a hole. This graph is 0(iV^^). 



It is worth to remind that in the theories with fermions the expansion of the fermionic 
effective action in Eq. ( [4. 2] ) is convergent only in the limit where the bare fermion mass is 
large, and therefore yields only limited information on the continuum limit of the lattice 
theory. It provides nonetheless a framework where analytical calculations can be performed 
and the planar equivalence can be tested explicitly. Note that for fermions in the adjoint 
representation, the full theory is still invariant under the center group Z^r and the large- A^ 
behaviour of the connected correlators is the same as the one we obtained above for the 
pure gauge theory. 
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Figure 2: Graphical representation of (tr {fl^))c.YM for = 3. The surface wraps N times around 
the thermal direction and ends in the gray area, which represents the integration of the product 
of N group elements, all oriented in the same direction. This integration yields a non-zero value, 
because the product of N fundamental representations of SU(A^) contains a singlet, given by the 
contraction with the completely skew-symmetric tensor. This graph yields a contribution 0{N'^). 




Figure 3: These graphs are the would-be leading contributions of respectively (tr ri)YM and 
((tr ri)^)c,YM- They are forbidden by the topology of the thermal dimension. 



4.1 Fundamental fermions 



When fermions in the fundamental representation are considered, the generic term of the 
r.h.s. of the Eq. (O) is: 



d"((P)p(pt)9; 



N'^-p-1 dc{oJi) ■ ■ ■ dc{ujn) 



Sf=0 



((P)P(Pt)g tr Ujcui) . ■ ■ tr ^K))e,YM 



(4.4) 



The expectation value that appears in the numerator of the r.h.s. of Eq. (4.4) is at most 
of order iV^-p-g-n^ since it contains p + q + n loops. Thus for n ^ 0, the whole term is 
at most of order N~^. Now, ii p = q the Yang-Mills contribution (n = 0) is 0(1) and 
dominates the sum. Otherwise, the leading contribution is at most of order A^~^: 



N'i-p-q 



)c,F 



O 



]S[2-p-q 
1 

iV 



0(1), \ip = q 



(4.5) 



For p ^ q, the 1/A^ scaling is only an upper limit for the asymptotic behaviour. 
However in the case of the expectation value of the Polyakov loop, the 0{N^^) is saturated 
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by the contribution: 



N ^ n\ ^ N'^-p-q dc(uji) ■ ■ ■ dc(uJn) 

n=0 LL}l...UJn 



Sf=0 



=^ E c(a;)(PtrC/(a;)),,YM + o(^) • (4.6) 

We recall that C(— 1) is the set of the closed paths on the lattice with winding number 
equal to —1 around the thermal dimension. 

From the large~A^ behaviour of the coefficients of the Taylor expansion, we see that, 
in the large- limit, the functional W is the same for Yang-Mills and for the theory with 
fundamental quarks: 

lim Wp{pL^a)= lim WVAf(«i'5) = 
N^oo N^oo 

^ lapP /Itr ni'^P) 

= l_y lim i^£iili£ . (4.7) 

^ (p!2P)2 Af-^oo iV2-2p V J 

This is the manifestation in this particular sector of the theory of the usual subleading 
contribution from the fermion determinant when fermions are in the fundamental repre- 
sentation. 

4.2 S/AS fermions 



Again, in the generic term of the r.h.s. of the Eq. (4.2): 

1 d'^{{py{p^r),^s/As 



]Sf2-p-q dc{uJl) ■ ■ ■ dc{uj„ 



(4. 

Sf=0 



each derivative in c{uj) produces the insertion of a vertex tr R[U{uj)]. When the fermions 
are in the (anti)symmetric representation, we can use the algebraic relationships: 

S/AS[!/] = <'^^'^f^'^^> . (4.9) 

In general, the insertion of some term of the form (tr [7)^/2 will disconnect the expec- 
tation value ((P)P(Pt)9) 

c,s/AS- In order to gain some confidence with these computations, 
we report in the following subsections the explicit computation of the leading orders of 
((-P)^)c,s/AS5 {PP^)c,s/AS ^'^d {P)c,s/AS- The end of the section is devoted to a discussion 
of the general case. 

Before entering into the details of the computations, we summarise here the results: 
j^^^z^^ = 0(1) , lip-q even , 

{{pnp^y)c,s/As ^ ( 1 



O -— , \ip — q odd and odd , 



, \ip — q odd and A^ even . (4-10) 



m'{p^r)c.^iA^ 



N'^-p-q 
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Hence the large- A'^ limit for the functional W yields: 
lim W5/^s(a,a) = 



^ pig! Af->oo N^-p-Q ^ ' 

p,q=0 
p—q even 



Since all the terms with even {p — q) contribute to the sum, we get 

Ws/As{a,a) = Ws/As{ua,ua) 

in the planar limit if and only if n = ±1. Therefore a Z2 symmetry is recovered in the 
large- limit. 

4.2.1 Computation of {P)s/AS 

Let us start with the single-derivative term in the expansion 



1 ^(J)s/AS ^1 R[Uiu;)])YM - (P)YM(tr R[U{co)])ym} = 

N dc{uj) 5^^o N 

=^{iP U{o,)]')ym ± {P tr [U{Lof ])ym} = 

=^ {(i^ [tr U{uj)f)c,YU + 2{P tr U (uj)) ^^Yuitr C/(l^))c,ym 

±(Ptr [U{Lof ])c,ym} . (4.12) 

Let w be the winding number around the thermal dimension of the closed path uj. We shall 
analyse each term in turn. The term {P [tr U [uj)]'^) c,ym saturates its behaviour only 
if 1 + 2^1; = 0, corresponding to correlators that are invariant under the action of U(l). But 
this equation has no integer solution, therefore the highest order is never saturated. The 
only non-vanishing contributions come from loops that satisfy 1 + 2w = kN (both k and 
A'^ must be odd). In this case, the first term goes like A^~^. 

The second term {P tr C/(a;))c,YM(tr U{u}))c,ym saturates its asymptotic behaviour 
only if each factor is separately invariant under the action of U(l). This corresponds to 
the equations 1 + w = and w = 0. Again, this equation admits no solution and the 
leading order is not saturated. We can look for subleading contributions. Since we want 
(tr C/(a;))c,YM not vanishing, it must he w = kN. Requiring that {P tr t/(u;))c,YM is not 
vanishing, we get 1 + w = 1 + kN = k'N which implies (k' — k)N = 1. This equation has 
no solution, and therefore the second term is always zero. 

The condition for the last term {P tr [U {uj)"^]) c^ym to saturate its A^'' behaviour is 
again 1 + 2w = 0. Looking for subleading contributions, we have to request 1 + 2w = kN. 
This equation admits a solution if both k and A^ are odd, as for the first term. But unlike 
the first term, the last one goes like A^~-^. 

Putting all these contributions together, the single-derivative in Eq. ( [4.2| ) is of order 
N~'^ and this behaviour is summarised by the formula: 



1 d{P)s/AS 



N dc{uj) 



± ^(P tr [C/(a;)2]),,YM + ^ ) • (4.13) 
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Consider now a generic term in the expansion (4^). We will not write the explicit 
expression, but it should be clear that it can be written as a sum of products of connected 
expectation values in the pure Yang-Mills theory. Let Wi be the winding number of the 
path uJi. Only products that are invariant under the action of U(l) saturate the leading 
contribution; hence the sum of all the winding numbers must vanish: 



l + 2^u;i = 0. (4.14) 



The factor 2 is due to the two-index representation. For each closed path uJi with winding 
number Wi, a loop tr R[U{lu)] with A^-ality equal to 2wi is inserted. The equation above 
has no solution, therefore the leading order for {P)c,s/AS is not saturated. This argument 
lets us to conclude that: 

li- = > (4-15) 

but we still cannot say if it is a 0{1/N) or rather a 0{1/N'^) as the single-derivative term 
suggests. Since this is not essential to our discussion of the center symmetry, we state 
simply the result as: 

!^=±ly]_y y y 

N N ^ n\ ^ ^ ^ 

n=0 k odd ■ZiCCi^ ^^^^ ) '^^ ■■'^^ 



c{ud)c{uji) ■ ■ ■ c{un) ^ y y ) \i<^ 



dc{uji) ■ ■ ■ dc{uJn) 
which means that {P)/N is of order 0{1/N'^). 



Sf=0 



+ (4.16) 



4.2.2 Computation of ((P)^)c,s/AS 

We want to check here that the leading term of ((-P)^)c,s/AS) which is expected to be A^*^, is 
actually different from zero. This can be easily shown by computing the single-derivative 
term in the expansion Eq. ( [4. 2D . Starting from the definition of the connected correlator: 
c,s/AS — (-P^)s/AS ~ (-f)s/AS' assuming that w is a generic closed path with winding 
number —1 around the thermal dimension, we can rewrite the single-derivative coefficient 
of the expansion of (-P^)c,s/AS ^-s: 



c,S/AS 



dc(Lij) 

_ 1 

~ 2' 



= {P^R[U ii0)])YM - {P^)YM{R[Uiu;)])YM 

Sf=0 



P2[tr U{u;)]^)ym 



= ^ {(P2[tr U{uj)f)c,YM + 2{P tr U{uj))c,ym{P tr U{uj))c,ym} = 

= {PtTU{u;)),^YM{PtTU{u;)),^YM + o(^-^^ , (4.17) 

as illustrated in Fig. ^. This is exactly a 0{N^) term. 
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Figure 4: This grapli represents a leading contribution of ((tr ^)'^)c,s/AS- The pair of blue loops 
represent the insertion of a vertex tr R[U{u!)], due to the derivative with respect to c{lu). 



4.2.3 Computation of (PPt) 



c,S/AS 



This is the easiest example: since PP^ is invariant under the action of U(l), the zero- 
derivative term (PP^)c,YM in the expansion ( [4.2| ) trivially saturates the A^*^ behaviour of 
(^'^'^)c,S/AS. 



4.2.4 The generic {PP{P^Y 



'c,S/AS 



As already discussed in the case of {P)c,s/ASj each term in the expansion ( |4.2| ) of the 
generic connected expectation value can be written as a sum of products of connected 
expectation values with respect the YM vacuum. Each of these connected expectation 
values contains some P's and P^'s from the original expectation value, and some tr U{uj) 
from the derivatives with respect to the coefficient c{uj). When the derivative with respect 
to c{LOi) is computed, a [tr [/(wj)]^ or a tr [[/(wj)^] is inserted. In both cases, if Wi is the 
winding number of the path uji, a parallel transport with winding number 2wi is inserted. 
A necessary condition for the leading order to be saturated is that the overall winding 
number, that is the sum of all the winding numbers of the operators involved, must be zero 
(this corresponds to the invariance under action of U(l)): 



P 



Wi 







(4.18) 



This argument implies that p — q must be even. If it is odd, the leading behaviour cannot 
be saturated. 

Now we want to show that this is also a sufficient condition, by explicitly constructing 
a term that saturates the leading behaviour. The p = q case is trivial since the leading 
behaviour is saturated by the pure YM term: 



c,YM 



+ ... 



(4.19) 



Consider now the p > q case (the opposite can be obtained by charge-conjugation). If u 
is a path with winding number —1, a leading contribution comes for instance from (see 
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1 d"((P)P(pt)'?)^ 




Figure 5: This graph is a leading contribution of ((tr Q)^tT ^^)c,s/AS- The pair of blue loops 
represent the insertion of a vertex tr R[U{uj)], due to the derivative with respect to c{lo). 

5. The effective potentials in the large— mass expansion 

The analysis of the center symmetry of the effective potential can be deduced entirely from 
the results presented in the previous section. Nonetheless, it is interesting to compare those 
results with the prediction in Ref. ||l^ about the effective potential of the Polyakov loop, 
bearing in mind that in Ref. |]l8| only the lowest-dimensional terms (relevant operators) 
are kept in the effective potential for the Polyakov loop T{z,z). They can be summarised 
as follows. 

Pure gauge or adjoint fermions. 

r(z, z) = a2zz + a4(zz)^ + a^iz^ + z^) (5-1) 

Fundamental fermions. 

r(z, z) =bi{z + z)+ a2zz + ai{zzf + aAr(z^ + z^) (5.2) 

S/AS fermions with N even. 

T{z, z) = 62(2^ + z'^) + a2zz + ai{zzf + a^iz^ + z^) (5.3) 

S/AS fermions with N odd. 

r(z, z) = hi{z + z) + 62(2^ + z"^) + a2zz + aA{zzf + aN{z^ + z^) (5.4) 
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In the large-A^ limit not all the coefficients survive. In particular: 



lim a7v(iV) = , (5.5) 

AT— »oo 



for all the theories; moreover, in the case of fundamental fermions: 



lim bi{N) = , (5.6) 

while in the case of S/AS fermions: 

lim bi{N) = , lim 62 (iV) ^ , (5.7) 

AT— KX) AT— >oo 

As shown in the previous Section, this behaviour is constrained by the symmetry- 
properties of each term in the potential. It is nonetheless intersting to check explicitly the 
last two equations. 

5.1 The 61 coefficient 

Since a{z,z) = —2dzT{z,z), from the Taylor expansion of the effective potential: 

6i = a,r(0,0) = -^ao , (5.8) 

where is the real source, defined by the condition z{ao,ao) = z{ao,ao) = 0, which is 
the same as: 

{da + da)W{ao,ao) = . (5.9) 
Using the series expansion of W{a,a), the above condition becomes: 

{trn) ^ (-«o)" ((Retr Q f+X 

^-"o>^(;^ = 0. (5.10) 

n=0 ^ ^ 

Since ^^j^ vanishes in the large- AT limit, the equation can be solved iteratively. At the 
leading order: 

an = ^ ^ — 7i f- subleadings . (5.11) 

N{{Retinf)c 

Expanding the denominator, the bi coefficient is: 

(tr ri) 

61 = -: r- + subleadings , (5-12) 

AT ((|tr + ((tr n)\j 

and it vanishes in the large- A^ limit, with both fundamental and S/AS fermions. 
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5.2 The 62 coefficient 



The second order in the Taylor expansion of the effective potential is obtained by inverting 
the Hessian matrix of the functional W: 



262 0.2 
a2 262 



(0) 



1 [ dlW dadaW 

'4 I d^d^w dlw 



The entries of the Hessian matrix in the large- limit are: 

1 ((tr 0)2 exp{-AraoRetr J]}) 



+ 



4 (exp{-iVaoRetr il}) 
1 /(tr exp{-A^QoRetr f]}) 



+ 
2 



1 



(exp{-A^aoRetr J]}) 



((tr 17)2), + ... 



(5.13) 



(5.14) 



1 (|tr exp {-A^gpRetr 0}) 
4 (exp {-TVaoRetr Vt}) 
(tr O exp{-iVaoRetr Vl}) 



+ 



(exp{-iVaoRetr 0}) 
(|tr + 



Computing the inverse of the Hessian yields: 

(|tr Jlp), 



02 
62^ 



[((tr 17) 



((tr 17)2 



|tr 17[2) 



2[((tr 17)2) 



2\ l2 



|tr 17|2),]2 ■ 



(5.15) 

(5.16) 
(5.17) 



In particular, the a2 coefficient is always of order 1; the 62 coefficient is of order ^ for 
fundamental fermions, while it is of order 1 for S/AS fermions. 



6. The effective potential in the deconfined phase 

We consider here the high-temperature regime. In other words, the space is taken to be 
X Si, where the extension L of the compact dimension is much smaller than Aq^^. 
Antiperiodic boundary conditions conditions are imposed for the fermion fields, so that 
the path integral can be interpreted as the partition function for a quantum field theory 
at finite temperature. 

If e™^ , . . . e*^^ are the eigenvalues of the parallel transport 17 around the compact 
dimension, the effective potential for v\^. . . vn can be computed in the one-loop approxi- 
mation. We refer to Ref. |jl^ for the details of the computation. At fixed A'^, the effective 
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potential for respectively OrientiQCD and AdjQCD are: 



Vor{v) 



1 



N 



N 



/(O, Vi - Vj) - 2Nf J2 fi"^^^ + + ^ 

i,j=l i<j=l 
N N 

XI /(^' '"i ~ '"j) ~ ^f'^ firnL, Vi - Vj + tt) 



(6.1) 



(6.2) 



where m is the mass of the fermions and the function / is defined in terms of the modified 
Bessel function K2: 



^ _ 1-^ 2 + {gx)'^K2{gx) cos{gv) 
J{x,v)-^2^ 4 



(6.3) 



9=1 



In the large-A^ limit, the function v{i/N) = Vi can be defined. It is related to the 
eigenvalue density p{v{x)) = [v'{x)]~^. Thus, in the large-iV limit, the effective potential 
is a functional of v{x): 



VOr{v) 
iV2 



1 



[/(O, v{x) - v{y)) - NffimL, v{x) + v{y) + tt)] dxdy , 



(6.4) 



^^^N^ ^JaJo ''^''^ ~ ""^^^^ ~ Nffi^L, v{x) - v{y) + tt)] dxdy . (6.5) 



In the large- AT limit, the center Zn acts on the function v{x) as: 

v{x) — >■ v{x) + a . 



(6.6) 



Clearly, the effective potential for the AdjQCD is invariant under the full center action. 
Instead, the effective potential for OrientiQCD transform as: 



VOr{v) 



1 

[/(O, v{x) - v{y)) - NffimL, vix) + viy) +2a + vr)] dxdy . (6.7) 



Since fix, v) has a 2tt period in v, the effective potential is invariant if and only if a = kir. 

This fact has consequences on the vacuum structure of the two theories. Since in both 
theories, the terms depending on vix) — viy) generate attraction between the eigenvalues, 
the minima are characterised by a vix) = v constant. In the AdjQCD v can take any value, 
instead in the OrientiQCD the fermionic term generates only two minima at = 0, tt, which 
correspond to P = ±1. 

Let us now discuss the implications for the orientifold planar equivalence. The charge 
conjugation symmetry acts on the gauge fields as: 



A^(a;o,x) -yl^(xo,x) . 



(6.8) 



The vacua of the OrientiQCD are invariant under charge conjugation. Instead, only two of 
the infinite vacua of AdjQCD are invariant under charge conjugation. These are precisely 
the two corresponding to P = ±1, and the orientifold planar equivalent must be valid 
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in this two vacua. It can be objected that, since all the vacua of AdjQCD are unitarily 
equivalent, the orientifold planar equivalence should be valid in all the vacua. However 
the matching between the observables in the two theories is more involved. Consider for 
instance a Wilson loop W^(r) = tr U{T)/N along the closed path T, wrapping w times 
around the compact dimension. Let {W(T))a^Ad} be its expectation value with respect to 
the vacua of AdjQCD identified hy v = a, and let (VK(r))or be its expectation value with 
respect to the vacua of OrientiQCD identified hy v = 0. The following equalities hold: 



(VF(r)),,Adj 



'(w^(r))o,Adj 



'{W{T))or 



(6.9) 



in the planar limit. Therefore the observable VF(r) in the a-vacuum of AdjQCD correspond 
to the observable e^°'^W{T) in the 0-vacuum of OrientiQCD. In the a-vacuum of AdjQCD 
an unbroken charge conjugation can be defined, by properly composing the naive one with 
the center symmetry. On the gauge field, it acts like: 



Ao(xo,x) 



.7-, , 2a 



Afc(xo,x) -Al{xo,x.) 
It acts on the Polyakov loop like: 



(6.10a) 
(6.10b) 



P(x) = tr Pexp 



e^*"tr Pexp 



I I Ao{xo,,:x.) dxo 


L 

t I Ao(a:o, ,x) dxo 




(6.11) 



and in particular it does not change the expectation value of the Polyakov loop. 



7. The center symmetry in the Hamiltonian formahsm 

It is instructive to reproduce the results of the previous sections using the Hamiltonian 
and coherent states formalism. This approach uncovers a deeper picture of the center 
symmetry, and its remnants in the large- limit. 

The large- A^ limit of gauge theories via coherent states was introduced in This 
formalism was used in ^] to prove orientifold planar equivalence, by defining an orientifold 
projection on the parent SO(A^) gauge theory with fermions in the adjoint representation. 
A review of these concepts is beyond the aims of this paper, therefore the reader should 
refer to the bibliography for details. 

The center symmetry is not a symmetry of the Hamiltonian. Indeed the center acts 
only on the temporal component of the gauge field, which is not a real degree of freedom 
of the theory: it is the Langrange multiplier for the Gauss constraint. Let us review some 
of the details here. 

If H is the Hamiltonian of the gauge theory on the lattice, the partition function is 
given by: 

Z = tr (e-^^P) , (7.1) 
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where P is the projector on the gauge-invariant states. If g is an element of the gauge 
group, let ^xb] denote the unitary operator, which acts on the Hilbert space by producing 
the gauge transformation g in the point x. As a consequence of the invariance of the Haar 
measure, the projector P and the partition function can be written as: 

p = / ll{g^[n{^)m^)} , (7.2) 

X 

^ = ytr (^e-^^llg^M^)]^ ]ldn{^) . (7.3) 
If {ipn} is a basis for the Hilbert space, the trace in the integral can be written as: 

tr (e-^^n^x[l^(x)]) = 5:(^„|e-^^|^r) , (7.4) 

\ X / n 

where V^i^^ is obtained by applying the gauge transformation 17 to the state ipn- From the 
equation above, r2(x) is the SU(A^) phase that the state ipn acquires after a translation 
around the temporal direction. By writing the matrix element of e~^^ as a functional 
integral, one sees that r2(x) is the parallel transport around the time direction and tr il(x) 
is the Polyakov loop. 

All the potentials of Sect. |5| can be written using the Hamiltonian formalism. For 
instance, the probability distribution of the parallel transport r2(x) is: 

e-snin)^^,^ |e-/^/^JJg^[0(x)]| , (7.5) 

If u is an element of the center, a center transformation is defined as Q{x.) uil(x). It 
does not affect the degrees of freedom in the Hamiltonian, but only the Gauss constraint: 



g-5n(nf2) 



itr \e-^^l[g^[un{^)]] . (7.6) 



When a center symmetry is present, it is not a symmetry for the quantum system in a 
proper sense: it is not implemented by a unitary operator on the Hilbert space. It is a 
symmetry only of the potential S^{Q). 

This fact implies that an analysis of the center symmetry in the large- limit using 
the coherent states formalism cannot be developed in a straightforward way. Indeed, in 
the coherent state formalism, the Gauss constraint is completely solved by taking only 
gauge-invariant observables as degrees of freedom. Thus, no analog of the Polyakov loop 
exists. 

This problem can be circumvented as follows. Rotate the system (in Euclidean space- 
time) and interpret the gauge theory on an infinite space and at finite temperature, as 
the same theory at zero temperature and on a space Si x with antiperiodic boundary 
conditions for the fermions. After the rotation, the Polyakov loop becomes the parallel 
transport around the spatial compact dimension, and the center acts now on the physical 
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degrees of freedom of the theory. If S is a plane orthogonal to the compact dimension, and 
is the set of all the positive links departing from sites of S and orthogonal to it, 
the center acts as: 

Ue^uUi iiieCti^). (7.7) 

Consider a Wilson loop Ty(r) in the spatial lattice. The center symmetry counts the 
winding number wiT) of the Wilson loop: 

W{T) ^ u'^^^^WiT) . (7.8) 

It is clear that, unless u (z Z2, the center mixes the real and imaginary parts of the Wilson 
loops, therefore it does not commute with the charge conjugation symmetry. 

Let us come now to the orientifold planar equivalence. The SO{N) parent theory has 
a Z2 symmetry in the large~A^ limit. This symmetry is mapped through the orientifold 
projection in the Z2 symmetry of the OrientiQCD in the large-A^ limit. Instead, the 
AdjQCD has a Z^r symmetry, but only the Z2 subgroup in the large- limit maps the 
neutral sector into itself. We conclude once more, that Z2 (and not Z^) is the only 
symmetry of OrientiQCD that we can deduce from the orientifold planar equivalence. 



8. Conclusions 



In this work we addressed the issue of the center symmetry of orientifold theories, by 
considering SU(A) gauge theories with fermions in the symmetric/antisymmetric two- 
index representations in the large- limit. Our approach is based on the idea that the 
right tool for an exhaustive analysis of the center symmetry is the quantum action (or 
equivalently the generator of the connected correlators) for the Polyakov loop in the theory 
at finite temperature. We investigated the quantum action in both the confined and the 
deconfined phases and we conclude that the orientifold theory is invariant under the Z2 
subgroup of the center. 

In the deconfined phase the invariance under (the spontaneously broken) Z2 implies the 
well-known existence of only two degenerate vacua. The analytic computation can be car- 
ried out in the high-temperature regime, where the one-loop approximation holds. In this 
case the effective potential of the Polyakov loop has two degenerate minima corresponding 
to {P)/N = ±1. 

In the confined phase the generator of the connected correlators of Polyakov loops can 
be analytically computed in the large-mass expansion of the lattice theory. We showed 
that all the terms, that are not invariant under Z2, vanish in the large-A limit. On the 
other hand the Z2-invariant terms are different from zero, showing that actually the center 
symmetry group is explicitly broken to the Z2 subgroup in the large- A limit. 

A previous analysis of the center symmetry of the orientifold theories can be found 
in Ref. |11|, where the authors concluded that the Z^- center symmetry is dynamically 
recovered in the large-A limit at zero temperature: the part of the fermionic action that 
is not invariant under the action of Z]\f decouples from the expectation values in the large- 
A limit in the confined phase. We want to point out that this result was obtained for 
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observablcs with trivial topology with respect to the temporal compact dimension. Only 
these observables are relevant at zero temperature. In particular the dynamics of Wilson 
loops is determined by the Zn invariant part of the fermionic action and this fact implies 
the stability of the A;-strings. However at non-zero temperature, the theory contains also 
observables with non-trivial topology with respect to the temporal compact dimension. 
We showed explicitly that the Polyakov loop couples to the part of the fermionic action 
with non-trivial winding number and this is in general a leading effect in the large-iV limit. 
Therefore, although the Z]\[ symmetry is dynamically recovered in a subsector, the whole 
orientifold theory is invariant only under the action of Z2. 

This picture fits well with the predictions of the orientifold planar equivalence. In 
the large-AT limit the orientifold theory is equivalent to an SU{N) theory with Majorana 
fermions in the adjoint representation in a neutral sector, defined by all the single-trace 
gauge-invariant C-even observables. The latter theory is invariant under the full Zn 
center symmetry. Although the two theories have different symmetry contents, we showed 
that this fact is not in contradiction with the orientifold planar equivalence. Indeed the 
equivalence is valid only between neutral sectors, thus we can expect the same symmetry 
in the two theories only after removing all the states and the observables outside of these 
neutral sectors. In other word, the orientifold planar equivalence implies the equality in 
the two theories of the symmetry subgroups that map the neutral sector into itself. In the 
case of the center symmetry, this subgroup is Z2 for both theories. 
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A. Lattice theory and large— N counting 



In this work we are interested in computing connected expectation values of products of 
closed loops in pure Yang-Mills theories. The Boltzmann weight e^"^ in the path integral 
can be expanded in a series in 1/A. The expansion of the Wilson gauge action produces 
a series of monomials of elementary plaquettes. The expectation values are computed 
by integrating each link variable over the group manifold, and non-vanishing results are 
obtained only when the plaquettes from the expansion of the action produce a tiling of 
a surface whose boundary is given by the closed loops. Each plaquette in this expansion 
appears with a factor of N. The rules for the SU(A^) group integration arc known in detail, 
and the result of integrating over the group manifold a generic product of matrix elements 
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can be found e.g. in Refs. |20, The integration contributes further factors of 

j dUUi,j,ul^^ = , (A.l) 



1 



/ 



]\f(^]\J"i — 1) ^^^i''^'^^hM^i2,m2^ji,h ~^ ^il,m2^j2,h^i2,mi^ji,l2) ) 

(A.2) 

dUUi^j^ ■ ■ ■ Uij^jj^ = ■ (A. 3) 



Each trace over colour indices (i,j,l,m) contributes a factor of A^. Collecting all con- 
tributions one obtains that each diagram is proportional to A^-^, where x is the Euler 



characteristics of the surface spanned by the tiling, see e.g. Ref. |16]. 
For the diagrams considered in this work: 

X = 2-2H-2B , (A.4) 

where H is the number of handles, and B the number of boundaries of the surface. The 
number of boundaries is given by the number of closed loops. One can readily see that the 
planar limit is obtained by considering surfaces with H = 0. 

Let us now describe in detail the computations that appear in the derivation of the 
results in Sect. |5|. 

Consider the connected expectation values: 

n 

(JJtr C/(a;i))c,YM . (A.5) 
1=1 

We will prove that the leading A^^^" behaviour is saturated if loops with winding number 
different from zero only appear in pairs with opposite winding numbers. 

The easiest case is (tr U{u)i)tr: U {uj2))y M ,c-, with wi = —W2- In this case the surface 
connecting the two loops is a cylinder, wrapping \wi \ times around the thermal dimension 
(see Fig. 0). Notice that it is essential that the loops have opposite directions, since the 
cylinder cannot be twisted because of the topology of the torus. Since the surface is planar, 
this diagram saturates the leading behaviour of the connected expectation value: 

(tr [/(wi)tr U{uj2))ym,c = 0{N^) with wi = -W2 ■ (A. 6) 

Consider now the more complex case of (tr f7(u;i)tr U{uj2)tT C/(a;3)tr U{llI4))ym,c with 
wi = —W2 and w^, = —w^. The planar surface is built by connecting each pair of loops 
with a cylinder with an hole, and than gluing the two cylinders through the boundary of 
the holes (see Fig. |^). The Euler characteristic of this surface is x = ~2 therefore the 
corresponding diagram behaves like A~^. 

The above procedure can be iteratively generalised to the case of an arbitrary number 
of pairs of loops. Moreover an arbitrary number of loops with zero winding number can be 
included simply adding holes to the surface. The resulting surface is always planar. 
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Consider again the generic connected expectation value in Eq. ( [A.5| ). If some loop oji 
with winding number Wi different from zero is not paired to a loop with winding number 
—Wi then the leading behaviour cannot be saturated. Since this fact is not crucial for our 
work, we only illustrate it in some simple cases. 

For n = 3 and wi = -0)2 = 1, W3 = —2, the connected expectation value 

(tr U{uJi)tr U{uJ2)tr U{uj^))ym,c 

has zero net winding number. However it is not possible to build a surface connecting this 
three loops because of the topology of the torus. One can ask if it is possible to use the 
four-link integration to get a leading contribution. We shall proceed in steps to prove that 
this is not the case. 

First of all, we can deform each loop by taking a contiguous plaquette from the ex- 
pansion of the Bolzmann weight and integrating the common link, as depicted in Fig. |6|. If 
tr {AU) is schematically the loop and A^A^^tr {U^B) is the relevant term of the action: 

NX-^ j dUii {AU)iT {U^B) = A^Hr {AB) . (A.7) 

In this way, the loop is replaced by the deformed one times a factor 1/A. We can use 
iteratively this procedure for each loop until the glued plaquettes shape three cylinders 
that end in three loops in the same spatial point (Fig. |^. Now we start to integrate one 
of the links of these three loops, by using the four-link formula in Eq. ( |A.2D . If we want 
to integrate the link [/, we can write the three loops schematically as tr {AU), tr {AU), 
tr {U'^ A'^U'^ A'^) (remember that they all are coincident). The result of the integration is: 

dUti {AU)ii {AU)ti {U^A^U^A^) = 
2 



tr {AAUA^) - _ [tr {AA^)f = 

. (A.^ 



A2 - 1 ^ ^ A(A2 - 1) 

2N 2N 



A2 - 1 A(A2 - 1) 

In the last line we used the fact that A is the product of all the links around the loop but 
U, and therefore it is a unitary matrix. We conclude that this kind of diagrams does not 
contribute to the connected expectation value. 

The last non-trivial case we want to illustrate is n = 1 and wi = N. In this case the 
connected expectation value (tr U{loi))ym,c is invariant under the center Zn, but not U(l). 
Therefore a non-zero contribution can be constructed only using the A-link integration 



formula for SU(A) in Eq. ( A. 3 ). By using the same construction as above, the loop tr U{uJi) 
can be deformed (without introducing extra N factors) into a loop wrapping straight in the 
thermal direction (see Fig. ||), that can be written schematically as tr [(AC/)^]. Integrating 
the hnk U: 

1 

det^ = l. (A.9) 



/ 



dUti [{UA)^] — ^j, ill. .. ipf^hji ^iNjN'^Ji---jN 
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Figure 6: The Eq. (A_/7) is represented. In the left side, the big loop tr (AU) and the plaquette 
tr (WB): in the right side, the deformed loop tr (AB). 



HnH ) 

^ ^1 i ^ 

Figure 7: The connected expectation value (tr U{u!i)tr U{u!2)tr U{u)^))ym,c with wi = W2 — 1, 
i/73 = —2. The four-link integration in the gray area makes vanish this contribution, as explained 
in Eq. (^). 

This is a subleading contribution since the leading term of the expectation value of a single 
loop is expected to be proportional to N . 
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